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Abstract 

In this paper we follow the lines of recent works to investigate systems of two 
coupled real scalar fields defined by potentials that describe periodic inter- 
actions between the scalar fields. We work with polymeric chains containing 
periodic interactions between the coupled fields, and we investigate the topo- 
logical sectors to obtain explicit soliton solutions and their corresponding 
energy. In particular, we offer an example that considers deoxyribonucleic 
acid (DNA) as a system of coupled fields, and we present the main steps to 
describe DNA as a polymeric chain belonging to the class of systems of two 
coupled real scalar fields. 
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I. INTRODUCTION 



This paper deals with sohtons in nonhnear systems of two coupled real scalar fields, 
described by periodic interactions between them in bidimensional spacetime. This is im- 
portant in nonlinear science in general, since there are many systems which are described 
by periodic interactions among their relevant degrees of freedom. The present work is a 
continuation of former papers [0-0], in which we have set the basic ideas we are now going 
to enlarge to include some specific systems that are described by periodic interactions. Our 
aim here is to investigate the presence of soliton solutions in systems belonging to a general 
class of systems of coupled real scalar fields, and containing periodic interactions. Evidently, 
we believe that such systems can engender interesting physical contents, and can be used 
to map relevant degrees of freedom in quasi-one-dimensional periodic chains usually consid- 
ered to model actual systems in condensed matter, organic chemistry and biological science. 
Within this context, in this work we use the term polymeric chains to name these systems 
generically. 

For simplicity, we consider the case of two fields, namely and x- Systems of this kind 
are usually described by a Lagrangian density that contains a potential U = U{(f),x), in 
general a nonlinear function of the real scalar fields. It is this nonlinearity that enlarges 
the scope of the problem, since it can be mapped to many interesting systems in nonlinear 
science. In this paper, in particular, we shall investigate periodic interactions that can be of 
great interest in applications to nonlinear science due to their topological properties, which 
we are also going to explore below. In this case, as an illustration we shall consider the 
deoxyribonucleic acid (DNA) as a system of coupled fields. In this first atempt to describe 
DNA as a system of coupled fields belonging to the class already introduced in 0, the main 
motivation is to shed some light on the way this can be done. Since the basic motivation 
of the present work is to introduce generic polymeric chains and to deal mainly with their 
mathematical intricacies, we postpone to a future work a more specific investigation of DNA 
as a system of two coupled real scalar fields. Here we recall another recent work |Q, in which 
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an alternate route to solitons in hydrogen bonded chains is introduced. 

The paper is organized as follows. In the next Sec. ^ we briefly review the method of 
searching for soliton solutions we have already introduced in There we introduce a 

general Lagrangian density and set the conditions the potential has to obey, such that the 
static solution of the second-order equations of motion are obtained by solving first-order 



differential equations. In Sec. |Tl| we examine examples of periodic systems. These examples 



illustrate how to search for soliton solutions in systems described by periodic interactions 



between two coupled scalar fields. In Sec. [I^ we consider DNA as a two-field system, and we 
show how to describe it via the class of system here considered. In this case, we start from 
the model proposed by Watson and Crick, in which DNA is formed by two polynucleotide 
strands that constitute the double helix held together by hydrogen bonds. We end the paper 
in Sec. 0, and there we present some comments and conclusions. 

II. A GENERAL CLASS OF SYSTEMS 

A general Lagrangian density describing a relativistic system of two coupled real scalar 
fields in bidimensional spacetime is given by 

C = U^<pd^<P + U^xd'-x - U{4>, x) , (1) 

where U = U{(f), x) is the potential, which specifies the particular system one is interested in. 
Our notation is usual: we are using natural units, in which h = c = 1, and the metric tensor 
5f"^ is diagonal, with g^^ = —g^^ = 1. In this case the fields are dimensionless, and time {t) 
and space coordinate [x) have dimension inverse of energy. The potential will be specified 
with two kinds of parameters: (a, 6, ...) as real and positive dimensionless parameters, and 
(A,/i, ...) as real parameters having dimension of energy. 

The above system leads to the following set of equations of motion 

dl^ _ 

and 



^ _ ^ + ^ = (3) 

dt"^ dx^ dx 

In the standard way ^ of searching for sohton solutions one consider static field configura- 
tions, and so (f) = (f){x) and x = x{^)- this case, the equations of motion become 

_ dU 
dx^ d(j) ' 

and 

d\ dU 



(4) 



(5) 



H^H^^ + H^H^^, (7) 



H^H^^ + H^H^x^ (8) 



As it was already stressed in let us now consider potentials that can be written as 

,,,, , 1 (ohV 1 (dnV 

where the functions H = H{(j), x) is a smooth but otherwise arbitrary functions of the fields 
(f) and X- Here the equations of motion describing the static field configurations become 

dx'^ 

and 

fx 

dx"^ 

where we are using if^ = dH/dcp etc. 

At first glance, it seems that we rewrote equations (||) and (H) in a more complicated 
form, but this is not so, as it was already shown In this case the energy of the system 

can be cast to the form 

/oo dfJ 
^dx — = if(0(oo), x(oo)) - if(0(-cx)), x(-oo)), (9) 

and this is the minimum value for the energy, which is achieved when we impose the condi- 
tions 



and 




(11) 



As we can easely see, solutions of the first-order equations (|T0|) and ( pJ]) also satisfy the 



second-order equations of motion (0) and since H{(f), x) is smooth. Therefore, for the 
general system (|l|), when the potential has the specific form @, the second-order differential 
equations of motion and (||) can be replaced by the first-order differential equations (p!oD 



We omit details here, but it was already shown j^J^ that every soliton solution this class 
of systems can comprise is classically or linearly stable. Furthermore, in the above class of 
systems the energy corresponding to static solutions is bounded from below, and gets to its 
minimum value given by (H). The function H{(j),x) can also be used to define topological 
sectors, as shown in 0: The topological charge, which is conserved, can be cast to the form 



and in this case it is identified with the energy of the corresponding static field configurations. 
However, since static field configurations should go to vacuum states, asymptotically, to 
mantain the energy finite, we see that the topological charge is nothing but the difference 
between two values of H{(f),x), calculated at, say, (0, x) and which represent two 

neighbour vacuum states in the (0, x) plane. 

III. POLYMERIC CHAINS WITH PERIODIC INTERACTIONS 

The investigations introduced in the former Sect. || are done on general grounds, and 
can be generalized to the case of three or more fields straightforwardly. Here, however, our 
motivation is to show how the general procedure we have already introduced works when 
investigating specific systems. In former papers we have already introduced several 
examples, in which we presented explicit soliton solutions. There we have investigated 



and dni). 




(12) 
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systems identified by potentials containing polynomial interactions between two coupled 
fields. In this section we focus our attention on periodic interactions, and consider potentials 
described by the sine and cosine functions. 

As we have already seen, the basic function we have to take under consideration is 
H{(f), x), and this is the understanding we shall now follow to investigate some explicit ex- 
amples of systems that contain periodic interactions between the two scalar fields. The 
following polymeric chains are introduced with two main motivations: Firstly, to deal with 
mathematical intricacies that naturally appear in these systems and, secondly, as a prepa- 
ration for possible applications to actual systems. 

A. Polymeric Chain Number One 

As a first polymeric chain containing periodic interactions, let us choose x) in the 
form 

Hi{4>,x) = A* cos(0) -|- u cos(ax) + A cos(0) cos(ax). (13) 

Other forms similar to the above one can be obtained by just shifting ^ + 7r/2 and/or 
X x+'''"/2a- In the present case, however, the first-order equations obtained from x) 
are 

^ + sin(0) [At + A cos(ax)] = 0, (14) 



and 

dx 



-I- asin(ax)[i^ + Acos(0)] = 0. (15) 



After some algebraic manipulations, the potential can be cast to the following form: 

1, 



X) = g[A'(l + a^) + 2(/.2 + u^a^)] + 



+ oK(« -l)-2r]cos(2(/)) + 
o 

+ liX'^il - a^) - 2v^a^\ cos(2ax) + 
8 



+ -Az/a^cos(0) + -X^cos{ax) — 

1 1 

— cos(20) cos(ax) — ttAz/o^ cos(0) cos(2ax) — 

- \\^{\ + a^) cos(20) cos(2ax). (16) 



Here we notice that in the above general function ([T3|) the two first terms correspond to 
decoupled self-interacting scalar field. The self-interacting fields have different amplitudes, 
controlled by the parameters [i and and different periodicities, the difference being gov- 
erned by the parameter a. If we set A = in (|13D the two first-order equations (0) and (|15D 
become 

^ + /isin(0) = O, (17) 

and 

^ + z/a sin(ax) = 0. (18) 
These equations can be integrated to give 

0(x) = 2arctane-^^, (19) 

and 

X(x) = -arctane"'"'''^, (20) 



which correspond to sine- Gordon systems [Q. 

The third term in (|l^) represents interactions between the two scalar fields, which is 
controlled by the parameter A. The general system contains several parameters, and presents 
the discrete Z2 symmetry, and for a = 1 and fi = u this symmetry becomes the Z4 symmetry. 
This last case seems to be interesting, since it describes two identical sine-Gordon systems 
(represented by the (p and x fields) interacting with each other. Another interesting case 
is obtained by setting /i = z/ = in Eq. ([TB|), and this will be investigated in the next 
subsection. 



If we set a = 1 and ji — v the potential becomes 

1 

4 
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C^i(</',X) = ^[A' + V]- 



^/x2[cos(2</.) + cos(2x)] + 
+ ^ V[cos(0) + cos(x) - 

- ^A/x[cos(20) cos(x) + cos(0) cos(2x) - 

- 7A^cos(2(/))cos(2x). 



(21) 



Furthermore, the first-order equations get to the form 



dcf) 



— + sin((/)) [/X + A cos(x)] = 0, 



(22) 



and 



dx 



+ sin(x)[Ai + Acos(0)] = 0. 



(23) 



This system presents an infinity set of singular points. For |A| < the singular points 
are at ((n7r, m7r), and alternate between stable and unstable points. For |A| > |//| there are 
more singular points, and they are located at (arccos(— ///A), arccos(///A)). At the particular 
value |A| = one gets to the picture that the system presents only stable or unstable points 
at the center of square cells with the four sides being continuous lines of singularities, as 
depicted in Fig. [1]. The following picture then emerges: for |A| < (|A| > |/i|) the intra 
(inter) chain binding is stronger than the inter (intra) one; for |A| = we have intra and 
inter chain binding equally stronge. 




(24) 



Fig. [1]. A cell of singular points in (0, x) plane for i?2(0, x) 
when II — u — X and a — 1. 
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We focus our attention at the specific case = |A|. To identify topological sectors in 
this system, we now use the topological current introduced at the very end of the former 
Sect. |I|. Here we verify that there is only one topological sector, with topological charge 
Qt = 4|A|. This topological sector represents any of the infinity possibilities of connecting 
the central point of the square cell to a point at its border, by a straight line. See Fig. [1]. 
A particularly simple soliton solution is obtained by setting x = 2ri7r and A = /i. In this 
case we get, from the above first-order equations, 

^ + 2Asin((/.) = 0, (25) 

which is the equation one gets in the sine-Gordon system, and this was already solved. 



B. Polymeric Chain Number Two 

As a second periodic system, we choose H{(j), x) in the following form 

^^2(0,x) = A cos(0) cos(ax). (26) 

This system is a particular case of the former system, and we treat it separately because it will 
be used in the next section. Here we also note that other similar forms are obtained by just 
shifting — i> + 7r/2 and/or x ~^ X + 7r/2a. For instance, the function H2 = Asin(0) sin(ax) 
can be considered as an almost trivial case of the investigation we are now doing, using the 
above H2. 

In the present case, however, the first-order equations become 

^ + A sm{(j)) cos(ax) = 0, (27) 



and 



- — h Aacos(0) sin(ax) = 0. (28) 
ax 



The potential can be cast to the form 



9 



[/2(0, x) = g A'[(l + + (a' - 1) cos(20) + 

+ (1 - a^) cos(2ax) - (1 + a^) cos(20) cos(2ax)]. (29) 

Note that the system presents Z2 symmetry, which becomes the Z4 symmetry for a = 1. 

This system presents an infinity set of singular points in the (0, x) plane. The singular 
points can be identified by: For n, m = 0, ±1, ±2, points at (2n7r, 2rmi /a) and at ((2n + 
1)7?, (2m + 1)11/0) alternate between unstable and stable respectively, and the points ((2n + 
l)7r/2, (2m + l)7r/2a) are all saddle points. See Fig. [2]. 




Fig. [2]. A cell of singular points in (0, x) plane for iJi(0, x). 
Thin lines represent topological sectors. 

To identify topological sectors in this system, we now use the topological current intro- 
duced at the very end of the former Sect. |T|. Here we use (^61) to verify that there are two 
distinct topological sectors: The first one connects adjacent unstable-stable points, which 
we name the paralell sector, since it joins points by a straight line paralell to the or x axis, 
with the corresponding topological charge given by = 2|A|; the second sector connects 
adjacent saddle-stable or unstable-saddle points, which we name the transversal sector, in 
which we have = Here we notice that the topological charge (and so the energy) does 
not depend on the parameter a; also, since Qj^ = 2Qip, we see that energy considerations 
favor the presence of solitons at the transversal sector. 

To find explicit soliton solutions, let us first investigate the paralell sector. In this case 
one sets x = ^nn/a, and so the above set of first-order equations ( P?] ) and ( pSD reduces to 
the single equation 
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^ + (-irAsm(0) = 0. (31) 

This is the sine-Gordon equation that was aheady studied in the former subsection. 

To investigate the transversal sector, we consider the simplest case: a = 1. Here, we can 
use = X + 2mi, and the first-order equations (|27| ) and ( pH] ) reduce to 

^ + i(-irAsin(20) = O. (32) 

See the former subsection for more details. It is interesting to point out that the transversal 
sector is connected by a straight line if and only if a = 1. We have being unable to find any 
explicit analytical solution in this transversal sector for a ^ 1. 



IV. THE DEOXYRIBONUCLEIC ACID 

In this section we shall consider deoxyribonucleic acid (DNA) as a system of the form 
we have been treating in the former sections. Before doing this, however, we recall that 
the DNA model proposed by Watson and Crick is formed by two polynucleotide strands 
that form the double helix, which are held together by hydrogen bonds connecting the 
adenine-thymine (A-T, with two hydrogen bonds) and guanine-cytosine (G-C, with three 
hydrogen bonds) bases. Two-field mechanical models for DNA consider each one of the two 
field describing bases in each one of the two strands, with the two strands being formed 
by harmonic coupling responding to the stacking of the DNA chain. Since the intrastrand 
distance between bases is at least twice as larger as the corresponding interstrand distance, 
such a mechanical view appears to be interesting, although one still has to account for 
interations between interstrand pair of bases governed by hydrogen bonds. 

This mechanical picture for DNA is usual, and the complete model in general differs 
only on the way one decribes interactions between interstrand pair of bases, as interestingly 
considered in , and in references therein. Since we are here interested in finding soliton 
solutions, we shall now follow the line of reasoning presented in [0] to model hydrogen bonds 
in conjugated pair of bases. Before going on, let us recall that the above mechanical picture 
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of DNA neglects inhomogeneities due to the base sequence in each strand, and asymetries of 
the two strands. Within this context, although the continuum-limit approximation appears 
to be somehow severe for DNA 0, when one discards inhomogeneities in each strand and 
asymmetry of the two strands, we realize that the continuum-limit approximation apeears 
to be not so severe anymore. 

If we now examine the steps introduced in 0] to get to the equation of motion, we 
recognize that the harmonic interaction present in each one of the two strands give rise, in 
the continuum approximation, to the second-order derivative along the DNA chain. Because 
of this, in the potencial there considered, we are just left with several terms, but all of them 
related to interactions between conjugate pair of bases that connects the two strands. This 
is all we need to jump from this mechanical model to our system of two coupled fields. 

Here we consider the fields and x describing bases in each one of the two strands in 
DNA. Moreover, since we are dealing with relativistic systems, we already have space and 
time coordinates entering with second order in derivative in the equations of motion. Then, 
the function H = H{(f), x) have already introduced is to be thought of as generating 
the potencial to describe interations between conjugate pair of bases connecting the two 
strands, in a way such that if we set to zero those interactions, we should end up with two 
independent strands, none of them containing self-interactions. We then recognize that the 
polymeric chain number two, already investigated in the former section, is the field system 
that appears to be more adequate to map the DNA model introduced in Ref. [0]. 

To see this more explicitly, let us consider the function 



which represents the addition of two functions H{(f), x) belonging to the polymeric chain 
number two already investigated in the former section. In this case we get the first-order 
equations 




(33) 



- — h sin cos X ~ ^ cos sin x = 0, 
ax 



(34) 
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and 

- — h /icos0sinx — z/sin0cosx = 0. (35) 
ax 

This system has an infinity set of singular points, and this is illustrated in Fig. [3]. In this 
Fig. [3] the central point is stable (unstable), and the other four points are unstable (stable). 
Here we see that there is just one soliton sector, which connects the central point to one of 
the other four points. The energy corresponding to this single soliton sector is given by 

Em = |yu - z^|. (36) 



(37) 



Fig. [3]. A cell of singular points in (0, x) plane for H^{(j),x)- 
The thin line represent the topological sector. 

We deal with the first-order equations ( ^41) and (|35|) to write 

^^^^t^ = (-/x + z.)sin(0 + x), (38) 

which is solved to give 

+ X = 2 arctan 6^"'^+'')^ . (39) 

Furthermore, we can also write 

c^(0 - X) . 
dx 

and now we have 



-(/i + z/)sin(0-x), (40) 



-X = 2arctane-(^+'^)^. (41) 
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The pair of solutions are 

= arctan 6^"^+^)^ + arctan e-^^^"^'' , (42) 

and 

X = arctan e(-^+'^)^' - arctan e~^^+''^''. (43) 
It is interesting to see that in this case the potential has the form 

— - (^/i^ + z/^^ cos 20 cos2x — 

— /^i/sin20 sin2x- (44) 

On the other hand, in |]7J an interesting model to map DNA is considered. There the 
author uses the following potential 

Uz{(p, (p') = {—B — 2/5) COSV9 cosv?' + 

+ (— -B + /9) sin (y9 sin — A cos ip — \ cos 99', (45) 

where B is the parameter associated with the hydrogen bond energy, (3 is the parameter as- 
sociated with the dipole-dipole interaction energy, and A is the coupling constant associated 
with the dipole-induced-dipole interaction energy. These are the parameters that appear in 
the DNA model considered in Ref. 0. 

Here we see that if one takes, for instance. 



^i = V2\Ib + 2/3+ V3/3(25 + p) , (46) 

and 



V2\ B + 2(3- J?>(3{2B + (3), (47) 



our potential (0) exactly reproduces the one considered in in the case A = 0. This result 
is very interesting, and shows that when the dipole-induced-dipole interaction is absent all 
the soliton solutions present minimum energy and are linearly or classically stable. 
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We now recall that the main motivation presented in for studying soliton solutions in 
DNA is to perhaps find mechanisms to explain duplication in DNA. Within this context, if 
one takes the reasonable point of view that stable solitons will hardly lead to mechanisms 
to explain duplication in DNA, we can very naturally improve the model by just adding the 
dipole-induced-dipole interaction term to the potential. In this case we see that this new 
potential is not of the form required to belong to the class of systems already introduced 
in Sec. |T|, and so we may perhaps open the possibility of introducing instability for guiding 
us toward mechanisms to explain duplication in DNA. To work with this more realist case, 
we can proceed as follows: We add the dipole-induced-dipole interaction term as an extra. 



infinitesimal contribution to the potential; then, we follow the procedure introduced in |T0 
to calculate corrections to the energy of the soliton solutions even when we do not know the 
solutions explicitly. Here we recall that the above issue was also considered in Ref. 0, and 
there the author introduced another approach, alternative. On the other hand, when the 
dipole-induced-dipole interaction term is turned off we see that the energy Em = \fJ^ — ^\ 
can be written as Em = V^(A+ — A_), with A± given by 



A± = p + 4d±2^36{2 + 6) , (48) 

where we have set 6 = f3/B. Furthermore, as we have already shown, in this case there is 
only one topological sector, and the corresponding energy is given by the above expression. 
Therefore, this value for the energy may be seen as a reference value, and can be used as 
a bound for energy considerations in investigating duplication in DNA. Furthermore, we 
could also consider the polymeric chain number one, and in this case we would introduce 
interactions in each one of the two strands, and these self-interactions would certainly im- 
prove the model considered in |0 to describe the DNA chain. Evidently, these issues are 
out of the scope of the present paper, and so we postpone to a future work the very specific 
investigation concerning the DNA polymeric chain. 
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V. COMMENTS AND CONCLUSIONS 



In this paper we have investigated a general class of systems of coupled real scalar fields. 
This class of systems is defined by 



where H = H{(f), x) is a smooth but otherwise arbitrary function of the two fields (p and x- Iii 
this case, the second-order equations of motion corresponding to static field configurations 
are solved by the following first-order equations 

dx 



and 



dx ^' 



The above set of first-order differential equations can be seen as a dynamical system, 
and so we can take advantage of all the mathematical tools available to dynamical systems 
to deal with it. We have also shown that static field configurations in this class of systems 
have minimum energy and are classically or linearly stable, and so the presence of soliton 
solutions will certainly play some important role in understanding physical properties of the 
system. 

Our investigation is done on general grounds, but we have introduced some examples 



of polymeric chains in Sec. |T|. There we have investigated specific systems, described by 
periodic interactions, with the motivation of showing how to deal with issues that naturally 
appear in the general procedure. As we have shown explicitly, periodic interactions lead to 
potentials that can be made periodic too, and the topological behavior is ease to be identified 
and can be used to guide us toward finding soliton solutions. 

This fact seem to be important in nonlinear science, since in this case there are many 
systems which can be described by periodic interactions among their degrees of freedom. 
For instance, as it was shown explicitly, we have found a model that maps the model used 
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in 1^ to describe DNA as a two-field system when the dipole-induced-dipole interaction 
term is neglected. This result seems to be interesting since it informs that without the 
dipole-induced-dipole interaction term, the continuum version of the model considered in ||^ 
presents stable soliton solutions, and stable solitons can hardly be used to model open states 
in DNA. Within this context, the presence of the dipole-indiced-dipole interaction seems to 
be a small correction that induces instability and so gives rise to a mechanism that may 
perhaps lead to open states, thus explaining duplication in DNA. We shall return to this 
issue in a future work, in which we deal specifically with solitons in DNA, owing to explain 
duplication via the presence of open states in the double helix model of Watson and Crick. 

As an ending comment, we would like to add that we can perhaps find other applications 
for the models investigated in the present paper. For instance, we believe that the polymeric 
chains we have introduced in this work can be used to map other systems, in particular the 
model introduced in to describe topological solitons in polyethylene. These and other 
related issues are presently under consideration. 
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